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SUSY EVENT GENERATORS FOR LINEAR COLLIDERS 



Frank E. Paige 
Brookhaven National Laboratory, Upton, NY 11973, USA 

The status of event generators for SUSY processes at future Linear Colliders is 
briefly reviewed. 

1 Introduction 

Parton shower Monte Carlo programs, reviewed by Sjostrand in the previous talk0 
are important tools for relating perturbative QCD cross sections to experiment. 
They also make possible more realistic predictions of the potential for future ac- 
celerators to study new physics. SUSY at the TeV mass scale is a well motivated 
extension of Standard Model. A Monte Carlo approach is needed not only to incor- 
porate properly the effects of QCD but also to handle efficiently the many possible 
chains of decays. Rather detailed treatments of SUSY are now included in three 
parton shower Monte Carlo programs: 

• ISAJET 7.44@ which incorporates ISASUSY by Baer and Tata; 

• PYTHIA 6.1P which incorporates SPYTHIA by Mrenna; 

• HERWIG 6.1P which takes part of its treatment of SUSY from ISAJET but 
adds especially i?-parity violation. 

This talk describes the physics assumptions made in these event generators. It also 
discusses a few recent developments, such as the incorporation of matrix elements 
in ISAJET. 

2 Minimal Supersymmetric Standard Model 

The starting point for all three SUSY generators is the Minimal Supersymmetric 
Standard Model (MSSM), that is, the SUSY model with the same gauge group as 
the Standard Model and the minimal particle content. Each chiral fL,R has a scalar 
partner fL.Ft.i and each gauge boson has a spin-1/2 gaugino partner. There are also 
two Higgs doublets with spin-1/2 Higgsino partners. 

The most general such model allows weak scale proton decay. The simplest 
solution is to impose symmetry under i?-parity, where R — (— l)3(s--f')+2S'^ jf jg 
conserved, then SUSY particles produced in pairs and the lightest SUSY particle 
(LSP) is stable. It is also possible to avoid proton decay by assuming only B or 
L conservation. TheJ!-violating couplings are Higgs-like, and some but not all of 
them must be smalllj If they are all small, then the only effect of R violation is 
to allow the LSP to decay, e.g., Xi — * QQQi Xi ~^ or Xi 99^- Larger 

couplings would also affect other SUSY decays; so far this has been incorporated 
only in HERWIG. 
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SUSY must of course be broken. It is possible to add gauge-invariant masses 
for all the SUSY particles by hand, but at the cost of 105 (or more) new param- 
eters. Presumably these masses should be generated by some sort of spontaneous 
symmetry breaking. This seems to be impossible using only the MSSM, so one 
must introduce a hidden sector to break SUSY at a mass-squared scale F and com- 
municate that breaking to the MSSM at a scale M. In supergravity (SUGRA) 
models, M - Mpianck = [SnGNy^^^ = 2.4 x 10^^ GeV so F - 10" GeV. In 
Gauge Mediated SUSY Breaking (GMSB) models, the communication is via Stan- 
dard Model gauge interactions, and M and F are much smaller. This can give a 
light {Mfj ^ 1 GeV) gravitino into which all other SUSY particles decay. 

Electroweak symmetry is broken through the Higgs mechanism, leaving five 
physical Higgs bosons, h^H^A^H"^. The gauginos and Higgsinos then mix to give 
four neutralinos and two charginos xf. In most cases, the X2: Xi are mainly 
gauginos, while the heavier states have Higgs-like couplings. The scalar fermions 
also mix, especially for the third generation. 

A random choice of the 105 MSSM parameters will produce flavor changing 
neutral currents and CP violation. These can be avoided by assuming that all 
phases and flavor mixings are absent. This produces what is called in IS A JET the 
MSSM parameter set: Mg, /x. Ma, tan/3, Afg,, M^^, M„^, Ml^, Me^, Mq^, Af&^, 
Mt^, Ml3, Mr^, At, Aft, Ar, Mq„ M,„, Mc^, Ml,, M^^ Mi, and Ma- This is a 
plausible starting point but has no good theoretical justification. 

More tractable parameter sets can be obtained by assuming some sort of uni- 
versality. If all scalar and all gaugino masses are taken to be degenerate at the 
GUT scale, then one obtains the minimal SUGRA model. Electroweak symmetry 
breaking is driven by the large top Yukawa coupling, leading to a four parameters: 
toq; ^1/2, ^Oj E^nd tan/3 = Vu/vd, plus sgn/i — ±1. This model is incorporated 
using a numerical solution to the renormalization group equations in ISAJET and 
HERWIG and an analytic approximation in PYTHIA. Recently optional additional 
parameters have been added to ISAJET so that one can study deviations from it. 

An alternative is the minimal GMSB model, in which the messenger sector 
contains 5 + 5 multiplets of SU{5). There are again just four parameters, the 
messenger scale M, A — F/M, N5, and tan/?. Gaugino masses are oc aAN^, while 
scalar squared masses are cx a^A^A^s, where A — F/M. The gravitino mass, or 
rather the lifetime for gravitino decays, is also a parameter. ISAJET also incorpo- 
rates a few additional GMSB parameters. 



3 SUGRA with tan /3 > 1 

If tan/3 is much larger than one, then effects such as — and — f^' mixing 
become important. This has been properly included in ISAJET since version 7.32. 
The most important effect is to increase the splitting of the third generation masses. 
In particular, the fi can become light, so that the only two-body gaugino decay 
modes may be X2 ~^ t^t^ and xt ~^ f^Vr- Th^e branching ratios may therefore 
be greater than 99%, as can be seen in Figure ffl _ 

A limited amount of work has been done on such cases for the LHCQ While 
measurements are possible, they are clearly quite difficult. Measurements at a linear 
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Figure 1: Gaugino branching ratios vs. tan/3 in the minimal SUGRA model. 
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Figure 2: Dilepton mass distribution with (solid) and without (dotted) matrix element. 



collider may be easier, at least relatively. E.g., one might measure masses using a 
threshold scan rather than from t's hadronic decays. This deserves more study. 

4 Matrix Elements in SUSY Decays 

Three-body decays are often important for SUSY. While matrix elements have 
always been included in the calculation of branching ratios, they have mostly been 
ignored in the generation of events. 

The newest version of ISAJET incorporates matrix elements for decays of the 
form A — > Bff. One does not want to duplicate all the code used to calculate 
branching ratios inside the event loop. However, the most general matrix element 
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Figure 3: Muon energy distribution with (solid) and without (dashed) bremsstrahlung and 

beamsstrahlung. 



has only spin-0 poles in the Bf and Bf channels and spin-0 and spin-1 poles in 
the // channel. The strategy is to save the relevant masses and couplings when 
calculating the branching ratios and then to reconstruct the matrix element with a 
single generic routine. The same strategy should work for HERWIG and PYTHIA. 

The effect can be quite significant, as seen from the dilepton mass distributiori 
from xfx2 ~^ ^^^Xi^^Xi ^tt the Tevatron shown in Figure ^. Nogiri and Yamadfl 
have also emphasised the importance of including matrix elements. Incorporating 
spin correlations in a general way is more difficult and has not been attempted. 



5 ISAJET Bremsstrahlung/Beamsstrahlung 

Bremsstrahlung has been included in ISAJET 7.44 using the Fadin-Kuraev dis- 
tribution function. Beamsstrahlung has also been implemented using the formalism 
of Chen, et alia The user must specify T and cr^. Both have significant effects on 
visible endpoint distributions, e.g., for e+e^ — > as can be seen from Figure^. 

This work is supported in part by the U.S. Department of Energy under contract 
DE-AC02-98CH10886. 
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